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Abstract 



The zero-temperature TAP equations for the spin-1 Ghatak-Sherrington model are 
investigated. The spin-glass energy density (ground state) is determined as a function 
of the anisotropy crystal field D for a large number of spins. This allows us to locate 
a first-order transition between the spin-glass and paramagnetic phases within a good 
accuracy. The total number of solutions is also determined as a function of D. 



PACS Numbers: 05.70.-a, 64.60.-i, 75.10.-b 
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Recently different authors have investigated some infinite-range spin-glass models 
which display both continuous and first-order transitions lines [ffj-B). In most cases, the 
first-order transition line starts at a tri critical point and extends down to T = 0. Whithin 
the replica approach, one should in principle locate this line according to the program 
proposed by Parisi [H] . However, it is easy to see that this is a very hard task, and has not 
been achieved so far. The difficulty has its origin in the fact that at a first-order transition 
at least two phases with different symmetries coexist with the same free energy but dis- 
tinct order parameters. In most spin-glass models at least one of these order parameters 
must be obtained within an accurate numerical study of the full Parisi treatment. 

Another method to study disordered systems was introduced in order to avoid replicas 
and is widely known as TAP approach [O. It is also well known that this method presents 
another kind of numerical difficulty, since the TAP equations have an exponentially large 
number of solutions, most of them related to metastable states . Nevertheless, it is 



possible to solve the TAP equations for the Sherrigton-Kirkpatrick (SK) model ]nj and 
get some useful information on the nature of the spin-glass phase [0,0 • 



Following SK, Ghatak and Sherrington JIB] introduced a generalized model to the 
case of integer spin Si = 0, ±1, . . . , ±S and including a crystal-field term. This model is 
an example of the systems mentioned in the first paragraph above. For S = 1 it has a 
first-order transition line that for some time was the object of some controversies [171-19] . 
Although its location is presently known within the replica-symmetric solution, it was 
not yet determined from a full Parisi solution. Recently, Feldmann and Oppermann |20| 
showed that a fermionic Ising spin glass is equivalent to the spin-1 Ghatak-Sherrington 
model. These authors also considered a one-step replica-symmetry breaking (1RSB) in 
order to locate the first-order transition line. For T = they concluded that the transition 
is located at D m 0.881J, where D represents the crystal field and J 2 /N the variance of 
the random couplings, a result already known to one of us |2T| . This should be compared 
to the replica-symmetric result which gives D « 0.899 J [T9]. In fact it results in a tedious 
algebra but one can show that a second step in the replica-symmetry breaking (2RSB) 



procedure gives D m 0.880 J |2l |. Thus, one may suspect that a full Parisi treatment 
would give a value for D/J very close to 0.88. It is thus natural to look for an alternative 
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treatment in order to check these results. 

The TAP equations for the Ghatak-Sherrington model were already obtained a few 



years ago by Yokota ||22j| . This author did an extensive numerical study of the transition 
at a particular temperature, namely, T = 0.2J. He showed that the first-order transition 
should be located at D/J = 0.85 ± 0.05. However, he did not search for the zero- 
temperature transition which turns out to be simpler to analyse. On the other hand, 
the TAP equations for the analogous fermionic Ising spin glass were recently obtained 
by Rehker and Oppermann ||. Following |?D| it is easy to see that there is an exact 
mapping between the corresponding equations for both models. Rehker and Oppermann 
|| have numerically studied their equations for T = and found that the transition line 
is located at D « 0.8J, without any mention to error bars. They also obtained equations 
which determine the number of solutions for the corresponding TAP equations at any 
temperature, but were unable to solve them numerically even for T = 0. 

The purpose of the present note is twofold. Firstly, we show that carrying out a nu- 
merical study with some reasonable numbers of spins, the TAP equations for the Ghatak- 
Sherrington model can be used to determine the transition at T = to a good accuracy. 
Secondly, we show that total number of such solutions can be determined as a function 
of D (or, equivalently, // in the fermionic-glass case |5J). 

We consider the Hamiltonian 

H=-Y,J V S i S j + DY,$, (1) 

(ij) « 

where each spin S{ (i = 1, 2, • • • , N) can take the values —1, and 1 and the summations 
are over all distinct pairs The random exchange couplings have zero mean and 

variance J 2 /N. According to Eqs. (3-4) from Yokota |]22|| , the TAP equations for this 
system can be written as 

2 smh((3hi) 



exp(pAi) + 2cosh(phi 



_ 2 cosh (fife,) 

Pi ~ exp((3Ai) + 2cosh((3h t ) [ } 

where 

hi = Jijmj - prm ^ Jf^Pj - m]) , (3) 

3 3 



and 



^ = d-U j24(pj-rf) , (4) 

Z 3 



where mj and pj are thermal averages fo Si and 5*?, respectively. 
At T = the above equation simplify to 



rrii = sga(hi)Q(\hi\ - D) , (5) 



hi = J i 
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where is the Heaviside step function. The energy density (in units of J) is given by 

/ = ~J^jJ2 Jijmmj + ^jY,Pi > ( 6 ) 

where pt = mf at T = 0. 

In the paramagnetic phase all rrii are equal zero and so is its energy density, irrespective 
of D. The spin-glass solutions can be found numerically. We have used an iterative 
approach to search for such solutions by rewriting the magnetization equations (|^) as 
J7i» in +i = ^ 7 ({n r ii,n})- The criterion adopted for convergence was 

1 N 

T7 \ m i,n+l - m i,n\ < 10 6 . (7) 
^ V i=l 

In the present case this method works finely and allows us to obtain as many spin-glass 
solutions as we could. That will not be so if we were working in a non-zero temperature 



regime as happens in the SK model ||23|| . We were thus able to improve the results 
obtained earlier [0], analysing systems varying from a few spins up to iV = 1000 spins. 
We have also varied the number of realizations of random interactions, iV^. Within each 
realization, the number of samples, N$, was determined from distinct initial conditions 
obtained as spin-glass solutions of the TAP equations for D = 0. For a given sample, we 
determine the solution with lowest energy density / m j n at D = and keep only solutions 
with energy densities such that|/// min — 1| < .05. The solutions thus kept, typically less 
than 5% of Ns, are then used as initial conditions to upgrade the solutions for a new value 
of D. This method allows us to obtain the spin-glass energy density as a function of D 
for each surviving sample. Finally, we averaged over the number of surviving samples and 



over the number of realizations. Fig. 1 summarizes our findings for the energy density. 
The error bars for the energy density are comparable in size to the symbols used on that 
figure, which give us much confidence on our results. We have also verified that for small 
systems (up to N = 200), the present method reproduces the results presented in Ref. ||, 
but with strong fluctuations. Thus, we find that the first-order transition which occurs 
when the spin-glass energy density becomes zero is located at 

D/J = 0.858 ± 0.008 . (8) 

This result improves the one found previously |J. Nevertheless, it is also remarkably 
different from those obtained within the replica approach p^9|-pT|. We have no sound 
explanation for this discrepancy between two seemingly equivalent methods as TAP for- 
mulation and replica treatment. We hope that other methods such as Monte Carlo sim- 
ulation, or exact determination of the ground state for finite systems, could help us in 
determining the zero-temperature transition in a definite way. 

The total number of solutions to Eq. @ can also be computed using the methods 
introduced by De Dominicis et al. JllJ or Bray and Moore and can be shown to 
give the same results. The latter method will be used in this note. Let us rewrite the 
magnetization equations as 

rrii = (p(hi) , (9) 

where (p(hi) = sgn(/ij)©(|/tj| — D). We also need to introduce the Edwards-Anderson 
order parameter defined as 

1 N 

<?=T7E™'- ( 10 ) 

iv 1=1 

Hence the total number of solutions (N s ) is given by 

fl r+l f+oo 

(N s ) =N dq Hdrm Hdhi5(Nq-J2 

JO J — 1 ■ J~ oo 



mf) 



n 



8(m,i - (p(hi)) [5(hi - ^ J, 



ijTTlj 



\ j^i i 

where (O) means the average of O over the random bonds. Introducing integral repre- 
sentations for the delta functions involving q and hi in the above expression gives 



(N s ) = N 



-Hoc rj\ r+1 r+oo r+oo An,, 

ioc ZlTl J —X ■ J—oo • J—oo ■ v27T 



J] 6{rrii - ^ (hi)) exp I NXq - A ^ m 2 + i ^ - -qmUvI 



ex P ~'E 



;i2) 



(ii) 



Performing the average over the Gaussian bond distribution, the last factor in flT2| ) be- 
comes 



cxp 



J 2 



w) 



cxp 



J2 J2 / 



(13) 



where we have neglected terms that do not contribute in the thermodynamic limit. The 
final factor in ( |13"D is simplified using the identity 



J 2 



2N 



;e y< 



N 



rrii 



00 dt 2 
e 

—00 



rrii 



Assembling the results ([13]) and ([14]) into (0), we obtain 



(14) 



/jy _ jy3/2 /"^ /" +0 ° ^ /" +i00 ^ /" +1 lT^ m f +0 °Y\dh [ + °° TT ^ 



AT 



x exp 



2 *2 + NXq - A X; ™- - ^ £ y 2 + % J2(hi + Jtmi)vi 



(15) 



In the thermodynamic limit, the above expression is dominated by the saddle point of the 
integrand with respect to the variables t, A and q. Thus, we have 



(N s ) « exp (iV> 



T) 



(16) 



where (fix is the saddle point of 



■\t 2 + Ag + lnH 



(17) 



and 5 is given by 



[° *L e -*/*i + 2 /■ 
^/27^ i- 



B / J da; 



-(x-Jt) 2 /2q-\ 



-D/J V27T 



for D > 0. For D < we recover the result for the SK model which gives 4>t ~ 0.1992 



11,12.24]. Therefore, we only need to determine 4>t for positive values of D. This 



is achieved numerically solving the saddle point equations for t, A and q. The result 
is presented in figure 2. It is interesting to note that <pT has a smooth behavior as a 
function of D, increasing from 0.1992, reaching a maximum around D rs 0.550J and then 
decreasing continuously to zero at D « 1.225J. The maximum is attained when the spin- 
glass phase presents a large number of spins in the S = state, whereas the remaining 
spins occupying the S = ±1 states are still in conflict due to frustration and randomness. 
Thus, the spin-glass phase may become more complex for intermediate values of D. As 
this parameter increases still further, eventually more and more spins prefer to stay in 
the S = state and finally the paramagnetic phase becomes the unique stable phase. It 
is also interesting to note that in the region where the first-order transtion is expected 
to occur the number of solutions is almost the same as in the Sherrington-Kirkpatrick 
model. 
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Figure Captions 



Fig. 1. Spin-glass energy density / as a function of D / J for N = 300 (A), 400 (□), 500 (o) 
and 600 (o) spins. The first-order transition is located when the curves cross the horizontal 
zero axis, since the paramagnetic energy density is zero for any value of D. 

Fig. 2. The logarithm of the total number of TAP solutions per spin, (p? = N~ l In (N s ), 
as a function of D / J. 
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